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if  of  tht  three  velocity  components  were 
t&f  flow  in  a  curved  rectangular  duct  with  an 
P|d||||||^u9TOflle.  The  results  provide  the  needed  in 


''t£K\djftelaF  analytical  models  of  internal  flows  with 
Xhtof  velocities.  Zn  addition,  new  analytical  models 
predicting  three  dimensional  rotational 


te  carved  passages ,  --with  inlet  total  pressure  distortion. 

bOCUOSn  'the  computed  results  and  the  experimental 
these  solutions  model  the  important 
iffUsnelnnal  flow  field  characteristics. 


:v' 


££*fl 


j'M 


vv; 


,*£ 


l:-  l> 


PROFESSIONAL  PERSONNEL 


Professor 


*.  Tabakoff 
S.  Abdallah 
C.  AtooGhutous 
C.  Cheng 
N.  Meade 
J.  Cupito 


Professor 

Graduate  Research  Assistant 
Graduate  Research  Assistant 
Graduate  Research  Assistant 
Graduate  Research  Assistant 
Mechanic 


ADVANCED  DEGREES  AWARDED 


M.S.  Degree  % o  M.  Malak,  Thesis  Title: 

"An  Investigation  Of  The  Trailing  Edge  Condition 
In  The  Finite  Element  Solution  Of  Inviscid  Flow 
In  Cascade  And  Single  Airfoil”. 

Ph.D.  Degree  to  S.  Abdallah,  Dissertation  Title: 

.  "An  InviSdd  Solution  For  The  Secondary  Flow  in 
Curved  Passages". 


tf. 


RESEARCH  OBJECTIVES 


The  main  objective  of  this  research  work  was  to  investigate 
the  secondary  flow  phenomena  analytically  and  experimentally. 

The  experimental  work  was  conducted  to  obtain  measurements  for 
the  three  velocity  components  using  Laser  Doppler  Velocimeter 

(LDV)  in  a  curved  duct  with  inlet  shear  velocity  profile.  The  purpose 
of  the  analytical  work  was  to  develop  a  formulation  and  a 
numerical  procedure  for  the  solution  of  internal  three 
dimensional- rotational  flow  fields. 

ACCOMPLISHED  WORK 

Since  all  six  technical  papers,  which  were  generated  under 
AfOSR  sponsorship  are  attached  in  this  report,  only  the 
important  contributions  will  be  summarized  here. 

The  first  phase  of  the  analytical  work  was  aimed  at  develop¬ 
ing  a  numerical  computational  procedure  for  inviscid  incompressible — 
rotational  flow  using  a  marching  technique.  The  governing 
equations  for  the  through  flow  velocity  and  vorticity  components 
and  for  the.  streamlike  functions  in  the  cross-sectional  planes 
were  derived  from  the  conservation  of  mass  and  momentum.  The 
nusnrical  results  were  obtained  for  the  rotational  flow  field 
in  a  curved  rectangular  duct  with  constant  cross  sectional 
area  and  curvature,  and  compared  with  available  experimental  data. 

This  work  was  published  in  the  AIAA  Journal,  Vol.  19,  No.  8, 

1981,  pp.  993-999  (Appendix  1} . 


SUMMARY  OF  SIGNIFICAHT  RESULTS 

/* 

— -pTwo  formulations  wore  developed  for  "nodeling  inviscid  three-  3 
dimensional  rotational  flow  fields  in  curved  passages .  The  first, 
for  a  very  efficient  marching  solution  with  hyperbolic  equations 
governing  the  development  of  through  flow  velocity  and  vorticity 
profiles  along  the  duct.  In  the  second  formulation,  a  new  approach 
was  developed  to  satisfy  the  conservation  of  mass  in  three  dimen¬ 
sional  flows  by  defining  two  dimensional  streamlike  functions  on 
fixed  orthogonal  surfaces  in  the  flow  fielcL^  The  first  formulation 
leads  to  a  faster  numerical  solution  in  which  the  influence  of 
the  pressure  field  on  the  three  dimensional  flow  can  be  sensed 
upstream  only  through  the  curvilinear  coordinate  system.  The 
second  formulation  on  the  other  hand  is  more  complex  since  it 
models  the  elliptic  influence  of  the  three  dimensional  pressure 
field.  Computer  time  savings  are  realized  through  the  two 
dimensionality  of  the  equations  for  the  streamlike  functions  and 
their  Dirichlet  boundary  conditions .  Agreement  between  the 
computed  results  and  the  experimental  data  was  very  good  in  both 
cases . 

V*  The  streamlike  function  vorticity  formulation  was  also 
tested  for  its  ability  to  model  viscous  flow  effects  and  their 
elliptic  influence  in  a  square  duct.^The  viscous  elliptic 
solution  predicted  the  influence  of  the  duct  on  the  flow  field 
upstream  of  the  inlet.  In  addition,  the  computed  results  were 
in  very  good  agreement  with  the  experimentally  measured  through 
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urine  do  10.L  method,  with  the 
(i.  (17)  awd  (IS),  b  the  numerical 
h  Bqs.  (14),  (16),  and  (17b)  are 


se,rkusedtosalvethehypetbolfcEqs.(6)aad 
ai  conditions  Bqs.  (S)  and  (9).  A  forward- 
otiet  k  uaad  to  advance  the  solution  in  9 
the  derivatives  with  respect  to  r  and  <  are  ap- 
nfns  central  flnksdUlsnnce  formulas.  Two 
Ural  wbtione  aw*  bo  satisfied  by  the  solution 
of  options  (6).  (7),  and  (li).  The  first  integral 
sets  the  alobol  condition  for  the  conservation 


where  A  k  the  dua  crow  sectional  area  and  Qk  the  volute 
flux  rate.  Due  to  truncation  error,  the  numerical  solution  for 
v  may  not  satisfy  condition  Eq.  (19)  exactly;  consequently,  a 
uniform  correction  is  introduced  in  the  source  term 
mnosendag  the  right-hand  side  of  Eq.  (12)  to  satisfy  Eq.  (19). 


'»  ^  S  \ 
^  ^  N  \  \ 
'  >  \  \  \ 
'  '  '  \  \ 


%)  Tetsskr 


of  At 
_  ft Mwt 
MitfM.  I5providmftow 


15b 

taht  velocity 
t  of  the  present 
tot  umiyiwwi  wnn  un 
tin  Raft.  10, 11,  tad  16.  Tbt  experimentally 
tiaftivriociy  profRaof  Raf.  iSiaraprodacadinFIg. 


ft  mMnpar  sacoad.  it  cm  bt  wat  from  this  figure  that  the 
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at  the  centerline.  In 
la  uiaerrad  that  the  variation  in  the  r 


direction  it  not  too  rigniflnmt.  The  inlet  vaiodtjr  profile  of 
Fig.  2  waa  found  too  noiey  to  be  need  dfaecdy  A  the  a— erical 
solution.  Therefore,  the  experimental  data  is  shnulatad  with  a 
parabolic  variation  in  the  *  direction  to  ohtaia  the  numerical 
results.  Due  to  the  symmetry  of  the  inlet  profile  the  com- 
potatioos  are  carried  oat  only  in  the  leear  half  of  the  Aset. 
The  results  in  Ref.  10  ware  also  prsaentad  A  the  lower  hall; 
while  Refs.  11  and  16  presented  their  remits  and  comparison 
with  the  experimental  data  in  the  upper  half  of  the  dact.  The 
governing  equations  were  nondimsnsionaMmd  before  the 
numerical  solution,  so  that  aU  remits,  except  the  velocity 
contours,  are  presented  in  nondimsnsionai  term.  The  duct 
inner  radius  A,  and  the  maximum  flow  velocity  at  inlet 
were  need  in  the  normalisation.  The  numerical  computations 
were  carried  out  in  double  precision  on  the  AMDAHL  470. 
The  results  presented  here  were  generated  using  an 
11x11x45  grid. 

The  results  are  presented  in  the  form  of  velocity  and 
secondary  vottidty  contours,  and  vectors  showing  the 
magnitude  and  direction  of  the  secondary  velocities.  The 
velocity  contours  at  the  J0-,  60-  and  90-deg  turning  angles  are 
shown  in  Figs.  3,  5,  and  7,  respectively.  It  can  be  seen  from 
these  Figures  that  the  rotation  of  the  velocity  contours,  which 
were  parallel  and  horizontal  at  inlet,  is  very  significant  in  the 
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the  divergence  of  the  vortkity  vector  is  identically  equal  to 
aero.  This  can  be  expressed  in  terms  of  u,  v,  and  w  and  (  as 
follows: 
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The  Elliptic  Solution  of  the 
Secondary  Flow  Problem 

7Mr  paper  preeentt  the  elliptic  solution  of  the  inviscid  incompressible  secondary 
flow  in  curved pafssgrs.  The  three-dimensk^ei  flow  field  is  synthesized  between  3 
aeW  of  orthogonal  nonstream  surf  acts.  The  two-dimensional flow  field  on  tech  set 
of  surfacm  It  coaOdend to  be  molting  from  e  source /sink  distribution.  The 
distribution  end  strength  of  these  sources  ere  dependent  on  the  variation  in  the  flow 
properties  norms i  to  the  surfaces.  The  dependent  variables  in  this  formulation  an 
the  velocity  components,  the  total  pressure,  and  the  main  flow  vortkity  component. 
The  governing  equations  in  terms  of  these  dependent  variables  are  solved  on  each 
family  of  surfaces  using  the  streamlike  function  formulation.  A  new  mechanism  is 
jmphmttntad  to  anehange  information  between  the  solutions  on  the  three  family 
surfaces,  resulting  into  a  unique  solution.  Inaddition,  the  boundary  conditions  for 
the  resulting  systems  of  equations  are  carefully  chosen  to  insure  the  existence  and 
uniqueness  of  the  solution.  The  numerical  results  obtained  for  the  rotational  in- 
stead  flow  in  a  curved  duct  are  discussed  and  compared  with  the  available  or¬ 
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Wilkinson  [§)  to  obtain  Made-to-blade  solutions.  The  same 
method,  ha*  also  been  used  by  Novak  and  Hearshy  (9]  and  by 
Kateanh  (10,  11]  to  obtain  meridional  and  Made-to-blade 
solutions.  The  assumptions  and  approximations  used  for  the 
stream  surface  shape  and  stream  filament  thickness 
distribution  in  the  5,  and  S2  solution  using  this  method  are 
m«lnty  similar  to  those  discussed  previously  in  connection 
with  the  matrix  method.  — 

Several  problems  are  encountered  in  computing  a  syn¬ 
thesized  three-dimensional  turbomachine  flow  field  from  the 
solutions  on  the  S|  and  S2  surfaces.  These  problems  are 
related  to  the  exchanged  information  between  the  two 
solutions,  concerning  the  stream  surface  shape  and  the  stream 
filament  thickness.  Novak  and  Hearshy  [9]  pointed  out  that 
the  St  filament  thickness,  as  calculated  from  the  blade-to- 
btode  solution,  is  not  constant  upstream  and  downstream  of 
the  Made  row.  This  is  in  contradiction  with  the  requirement 
that  the  flow  must  be  treated  as  axisym metric  in  these  regions. 
They  also  discussed  the  effect  of  the  lean  of  the  S2  mean 
stream  surface  extending  upstream  and  downstream  of  the 
blades,  where  the  net  angular  momentum  changes  must  be 
zero.  Stuart  and  Hetberington  (12]  tried  to  use  a  technique 
similar  to  Wu’s  in  their  solution  for  rotational  flow  in  a  90- 
deg  bend,  by  synthesizing  the  three-dimensional  flow  field 
through  the  iterations  between  two-dimensional  solutions. 
They  were  unable  to  reach  convergence  in  the  iterative 
numerical  solution  and  had  to  abandon  this  approach.  They 
speculated  that  the  information  conveyed  between  the  two 
solutions  was  not  sufficient  to  produce  convergence. 

Th*  assumption  that  the  5,  stream  surface  is  a  body  of 
revolution  was  common  in  all  applications  of  Wu's  theory  in 
turbomachines  P-12).  This  assumption  is  valid  if  the  flow  is 
hrotational.  In  general,  the  5,  surface  is  twisted  and  skewed 
due  to  the  presence  of  the  secondary  velocities.  These 
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tram  vase  velocities  exist  due  to  the  streamwise  vortkity 
which  is  generated  by  the  turning  of  a  nonuniform  flow  with  a 
vortkity  component  in  the  curvature  plane  (13,  I4|.  The 
nonunifonnitks  of  inlet  velocities  in  turbomachines  are 
associated  with  the  hub  and  the  tip  boundary  layers.  Stream 
surface  warpage  pose  additional  mathematical  difficulties  in 
the  solution  of  the  rotational  flow.  Fagan  (13)  could  not 
obtain  a  solution  for  highly  rotational  flow  in  curved  ducts 
using  Wn’s  approach.  He  resolved  the  problems  encountered 
when  the  stream  surfaces  warpage  approaches  90  deg  through 
using  a  corkscrew  coordinate  system  that  operates  at  a 
specified  angular  rate. 

In  reference  1161,  Abdallah  and  Hamed  developed  an  ef¬ 
fective  method  for  the  solution  of  three-dimensional 
rotational  flow  in  curved  ducts,  in  which  the  through  flow 
velocity  and  vortkity  components  were  computed  using  a 
marching  technique  in  the  main  flow  direction.  This  led  to  a 
very  efficient  numerical  solution  whose  results  compared  a 
favorably  with  the  experimental  data  for  duct  flows. 
However,  because  of  the  marching  technique  used  in  com¬ 
puting  the  through  flow  velocity,  the  influence  of  the 
downstream  conditions  on  the  flow  characteristics  is  not 
modeled.  This  effect,  although  not  significant  in  duct  flow, 
may  be  quite  important  in  turbomachinery  applications. 
Barba  and  Langston  [17]  contrasted  the  Made  row  and  duct 
flow  problems  and  discussed  the  importance  of  the  eiliptk 
solution  to  the  flow  in  blade  rows. 

This  investigation  represents  an  eiliptk  solution  for  the 
internal  non  viscous  incompressible  rotational  flow  in  curved 
passages.  The  streamlike  function  method,  which  was 
developed  by  the  authors  [18]  for  the  efficient  numerical 
solution  of  the  continuity  and  rotationality  equations  is 
implemented  in  the  present  problem  formulation.  The 
dependent  variables  in  this  formulation  are  the  three 
stramlike  functions,  the  total  pressure,  and  the  throughflow 
vortkity  component.  The  equations  of  motion  are  satisfied  on 
three  arbitrary  sets  of  orthogonal  surfaces.  On  these  surfaces, 
two-dimensional  Poisson  equations  are  derived,  for  the 
streamlike  functions,  with  source  terms  representing  the  flow 
three-dimensiofiality.  The  source  terms  are  dependent  upon 
the  variation  of  the  flow  properties  in  the  direction  normal  to 
these  surfaces.  Because  of  the  dependency  of  the  solution  on 
each  set  of  surfaces  on  the  solutions  for  the  remaining  two  sets 
of  surfaces,  an  iterative  process  is  involved  in  the  solution. 
The  three  flow  velocity  components  are  determined  by  the 
source  terms  and  the  three  streamlike  function  derivatives. 
The  total  pressure  and  through  flow  vortkity  are  computed 
from  Bernoulli  and  Helmholtz  equations,  respectively. 

The  initial  and  boundary  conditions  for  a  closed  system  of 
equations  are  carefully  chosen  to  insure  the  existence  and 
uniqueness  of  the  solution.  The  no-flux  condition  at  the  solid 
boundaries  leads  to  Dirichlet  boundary  conditions  for  the 


streamlike  functions.  Downstream,  the  derivatives  of  the  flow 
properties  in  the  throughflow  direction  is  set  to  zero.  A 
Poisson  type  equation  with  Neumann  boundary  conditions  is 
derived  and  solved  for  the  static  pressure  at  the  inlet  plane, 
which  is  then  used  together  with  the  inkt  velocity  profile  to 
determine  the  inkt  total  pressure  profile. 

Numerical  resutls  are  obtained  for  the  case  of  rotational 
flow  in  a  curved  duct  with  rectangular  cross  sections.  The 
results  are  discussed  and  compared  with  the  experimental 
data. 


Matbeaaatkal  Formutmtioa 


For  simplicity  and  to  be  Me  to  compere  with  existing 
experimental  results  -in  ducts  [19J,  the  cylindrical  polar 
coordinates  are  used  in  the  following  presentation.  The  bask 
equations  for  non  viscous  incompressible  flow  are  expressed  in 
terms  of  the  three  velocity  components,  the  total  pressure  and 
the  throughflow  vortkity  components  as  follows: 

Conservation  of  mass 


1  d  ,  .  I  to  a*  . 
Conservation  of  momentum  in  /--direction 


U) 


r  99  1  *  dr 


Conservation  of  momentum  in  z -direct ion 


(2) 


*  r 1 


ill  ~ 
9z  J*  9z 


(3) 


where  P  is  the  total  pressure  divided  by  the  density,  and 
(tt.tt.w)  are  the  three  velocity  components  in  (rAzHiirections. 
The  throughflow  vortkity  component,  $,  can  be  written  in 
terms  of  the  cross  velocities,  u  and  w,  as  follows 


9u  9w 
~9z  ~  3r 


(4) 


BernoaUTs  Equation. 

Bernoulli’s  equation  is  used  instead  of  the  momentum 
equation  in  the  0-direction. 


u 


—  +  JL  — 

9r+  r  99 


-0. 


W 


Helmholtz  Equation. 

di  o  9(  91 

u  4.  _  4.  w 

9r  r  99  9z 


9v  {  9v  .  9v  I  ,  „ 

,,1T  +  710  + 


(6) 


.  Nomenclature 


duct  cross-sectional 


area 

fi,Zi  -  arbitrary  integration 

C  ■  contour  enclosing  A 

reference  points  on  the 

component  in  0- 

dC  «  incremental  distance 

r-  and  z -coordinates, 

direction 

along  C 

respectively 

(Ar,A$,Az)  * 

space  increments  in 

n  m  outward  normal  to  the 

S,  ,S2  m  blade-to-blade  and 

(r,0,z)  directions 

contour,  C 

meridional  stream 

P  a  total  pressure  divided 

surfaces,  respectively 

Subscripts 

by  the  density 

(u,u,w)  m  velocity  components 

h  * 

horizontal  surfaces 

P,  a  inkt  static  pressure 

in  (r,9j)  directions, 

/  a 

inlet  conditions 

divided  by  the  density 

respectively 

V  a 

vertical  surfaces 

m  cylindrical  polar 

X  *  streamlike  function 

C  » 

cross-sectional  sur¬ 

coordinates 

|  a  main  flow  vorticity 

faces 
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(13).  to  obtain  the  fottowing 
far  x*.  x*.  and  x».  The  integration  reft 
far  the  internals  fa  eqt£®Joae  (7-9)  « 

r,  m  R,  mif,  m  0,  respectively, 


Nrfl 
hm  to  be 


cm 

(26b) 

(260 

(TO 

<27b) 

(270 

(28«) 

(286) 

(28c) 

(290 

(296) 

(290 

(300 

(306) 

,(300 


* 


--=.-0 


(246) 


(24c) 


to  umO  to 

and*,: 

****t?M?-  « 

(21-25)  art  wad  together  with  aquations  (7-9) 


The  governing  equations  (5),  (6).  (10),  (1 1)  and  (12)  with  the 
conditions  given  by  equations  (14),  (13),  and  (26-30)  form  a 
dosed  system  which  is  solved  for  the  variables  x*.  x*.  X».  f. 
and/*. 


The  results  of  the  computations  of  the  secondary  flow  in 
curved  ducts  caused  by  total  pressure  inlet  distortion  are 
presented.  The  iterative  solution  procedure  is  based  on  the  use 
of  successive  over  relaxation  method  for  the  solution  of  the 
three  streamlike  function  equations,  and  Lax's  [20]  marching 
scheme  for  the  total  pressure  and  through  flow  vortidty 
equations.  The  results  are  presented  and  compared  with  the 
experimental  measurements  of  Joy  (19)  for  the  flow  in  a 
curved  duct  with  constant  curvature  and  rectangular  cross 
section.  Joy  (19)  obtained  flow  measurements  in  a  curved 
rectangular  duct  of 0.125  x  0.25-m(3  x  10-In.)  cross  section, 
0.375-ra  (13-in.)  mean  radius,  and  90  deg  turning  angle.  A 
large  velocity  gradient  was  produced  in  the  experiment  using 
screens  placed  before  the  curved  portion  of  the  duct,  which 
resulted  in  a  nearly  symmetrical  velocity  profile  at  inlet  to  the 
bend.  The  velocity  contours  in  the  lower  half  of  the  duct  are 
shown  in  Fig.  2.  The  computations  were  carried  out  in  the 
lower  half  of  the  duct  to  take  advantage  of  the  symmetry.  The 
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MNb  of  the  aMvatr:>w  are  compared  with  Joy’s  flow 
MMM  at  30,  40,  aod  90  dag  tontine  angles.  The 
qspdtiidasof  the  flcwrwigeMwfa  Ftas.2-3  are  nonulised 

dan  0mm  flgares  that  the  compwtad  results  are  in  good 
KIMM<  aih  i he  aspailmaatal  inmsiimiuivi  at  thr  TIT  anil 
99  d«  (anting  eagles.  There  is  a  lack  of  agreement  however 
*lgi  On  onperiateata!  remits  at  the  00  deg  tnrnhig  angle  inner 
waB  near  the  daet  centerline.  Other  investigators  (12,  13] 
aamparing  with  the  saase  eaperhaental  resuhs  speculated 
ehteaae  hoandary  layer  separation  there  and  attributed  the 
Wk  of  SgrssMsnt  to  h.  The  analysis  otherwise  prediett  the 
UMaeured  contour  rotation  caused  by  the  secondary  flow 

00  and  90  deg  Hinting  angles  are  presented  in  Figs.  0.  These 
contours  tiww  that  the  static  pressure  is  not  constant  over  the 
cross  sections  even  at  xaro  dm  turning  angle.  The  static 
pMMmgradtaat  hi  the  radkd  direction  is  necessary  to  balance 
the  centrifugal  forces.  Joy  (19]  did  not  obtain  static  pressure 
HMasurenMats  that  could  be  compared  with  the  present 
make;  however,  static  pressure  contours  similar  to  those  of 
Fig.  ie,  were  laeaeurnrt  by  Braun  [20]  in  a  curved  duct  at  zero 
taming  angles.  The  effect  of  the  secondary  flow  development 
on  the  pressure  distribution  is  demonstrated  in  these  figures 
by  the  variation  hi  the  shape  of  the  contours  with  the  duct 
turning  angle.  The  computations  were  carried  out  using  a 
aaHana  ft  x  9  x  31)  grid  hi  the  r-,  end  ^directions, 
respectively.  The  convergence  of  the  iterative  procedure  was 
very  fast,  as  shown  in  Fig.  7,  which  presents  the  average  of  the 
absolute  value  of  the  error  hi  the  calculation  of  «*,  the 
through  velocity  gomponint  from  the  x»  solution.  The  ex- 
change  of  Information  from  the  solutions  oa  the  vertical  and 
crom-sectioaaf  surfaces  to  this  surface  does  not  start  until  the 
seaand  iteration,  which  leads  to  maximum  error  that  con¬ 
sequently  decreases  very  rapidly.  The  solution  was  obtained 
using  double  precision  with  30  outer  iterations  and  120 
iMratiaas  hi  the  solution  of  the  differential  equations  for 
x»,  and  x,  on  el  49  surfaces  and  required  4.3  min  on  AM- 
ShjHL  3m  The  authors  have  not  attempted  to  optimise  the 
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CPU  time  through  changing  the  number  of  iterations  in  the 
solution  of  Poisson’s  equations  with  the  outer  iteration  cycles, 
or  to  reduce  the  CPU  time  through  the  use  of  noniterative 
methods  (21,  22]  in  the  code.  At  present,  direct  Poisson 
solvers  codes  are  being  studied  for  incorporation  into  the 
numerical  solution.  This  is  expected  to  lead  to  considerable 
CPU  time  saving!  when  it  is  combined  with  the  streamlike 
function  formulation  and  its  corresponding  Di  rich  let 
boundary  conditions. 


It  can  be  concluded  that  the  present  analysis  can  predict  the 
in  viscid  secondary  flow  development  caused  by  inlet  total 
pressure  distortion  and  the  results  of  the  computations 
compare  with  the  experimental  measurements.  Through  the 
elliptic  solution,  the  influence  of  the  downstream  conditions 
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mttuoMl  flora.  bancor  and 
hsvelogad  a  vary  aiaple  nodal  for 
'•inti  of  oligraptOc  flow  fialds  by 
velocity  vector  Into  irro-  ; 
and  rotational  parta.  they  uaad 
function  to  daaoriba  tba  rotation- 
for  tbia  special  olaaa  of  pro- 
l  Ahiflh  tba  rothalpy  and  antropy 
>  aim  orthogonal  to  tba  vortieity 


natbod  for  tba  alliptio  aolution  of  tbraa  • 
dinanaional  rotational  inconpraaaibla 
iaviacld  flow  in  which  tba  valocity  vac tor 
Ac  ayatbaalaad  from  two  dinanaional  aolu-  ! 
tiona  of  tbraa  atraanlika  funetiona  on 
tbraa  aata  of  orthogonal  aurfaoaa.  Aaida 
fron  inconyraaaibility,  no  additional 
■jpunyt  Itiar  warn  iaposod  that,  night  linit  ; 
tan  analysis  to  a  spaeial  olaaa  of  pro- 
blana.  Tba  prasant  work  rapraaanta  an 
aa tana Ion  of  that  analysis  to  cooprossible- 
flow  fialda  with  generalised  channel 
gaonetrias  using  orthogonal  curvilinaar 
body  fitted  coordinates.  The  problea  for- 
anlation  loads  to  a  vary  efficient  aolution 
since  tba  governing  aquations  for  the 
otBoCbllhc  funa tiona  consist  of  poiason'a 
equations  with  Dirichlat  boundary  condi- 
tiona  abila  tba  gowning  aquations  for 
total  pressure,  total  enthalpy  and  straan- 
WfOS  vortieity  are  hyperbolic.  Tba  con- 
vaegaaea  of  this  thraa-dinanaional  solution 
ia  vary  feat**  and  doss  not  suffer  fron 
the  pnahlnra  aaeountarad  with  the  tradi- 
•focil  qpaasi-threo-dinanaional  nathods  in 
tba  pcaaaaoa  of  straaawise  vortieity. 


coord inat as,  respectively,  p  is  tha  flow 
density  and  ,  hj,  hj  ara  tha  as  trie  coef- 

ficianta  of  tha  orthogonal  curvilinaar 
coordinates. 

Tha  aquation  of  conservation  of  aass 
la  identically  satisfied  through  tha 
definition  of  three  atraanlika  functions1-7 
Xjj,  Xy  *od  xc-  The  atraanlika  function  xh 

ia  defined  such  that  its  derivatives  on  the 
surface,  x-  -  constant  ara  related  to  tha 
two  velocity  coaponents  Vlh,  V,h  as  follows: 


:  "  hih3e  v2h 


*  l  fe  "*2  ,2*» 


I  ?3q  (h3xh)  “  “  Vlh  * 

The  following  definition  of  the  second 
atraanlika  function  xv  i*  given  in  terns 
of  its  darivativaa  on  the  surface  x,  -  cc 
stanti  * 


lx^  fh2xv*  “  “  hlh2pV3v 


(h.h.pV,.  )  dx,  (3a) 


IxT  ,nln3pv2h 


Analysis 


Tba  governing  aquations  uaad  ia  tha 
aaaarical  aolution  of  tba  oonpressibla 
three  dinanaional  i aviaaid  rotational 
internal  flow  ara  derived  fron  tha  basic 
aquations  of  conservation  of  aaaa  noawn tun 
and  energy  la  orthogonal  curvilinear  body 
fitted  coordinates.  In  the  analysis,  hyper* 
to lie  aquations  ara  derived  for  tha 

through  flaw  vortieity,  tha  total  enthalpy, 
while  elliptic  aquations  ara  derived  for 
three  seta  of  atraanlika  functions  which 
are  defined  on  fixed  orthogonal  surfaces. 
Tha  details  of  tha  analysis  for  inoonpres- 
aibla  flow  can  bo  found  in  reference  16. 

In  tha  following  derivations,  tha  sane 
approach  ia  applied  to  oonpressibla  flow 
using  orthogonal  curvilinaar  coordinates. 


rq  (h2xv>  *  h2h3pVlv  (3b> 

Tba  third  atraanlika  function  x_  i*  similar¬ 
ly  defined  aa  follows: 

:  fej  (hlxc>  *  ■  hlh3pV2c 
x2 

*  l  fsT  <hih2pV3v)  te2  (4a) 

X  •>  J 


-  »»il‘2pV3c 
x3 

*  [  hz  <hlh3pv2h)  dx3  (4b) 


The  aquation  of  conservation  of  aass 
for  ooaprassihla  flow  in  curvilinaar 
coordinates  is  given  by 

»2»,»V  ♦  fej  'WV 

.  Jij  <*,VV  -  0  (II 

where  V^,  Vj  and  Vj  ara  tha  valocity  con- 
ponanta  in  tha  directions  of  tha  x.,  x,f  * 


The  integrals  on  the  right  hand  sides  of 
aquations  (2)  through  (4)  represent  source 
terns  which  are  dependent  upon  the  flow 
crossing  tha  three  sate  of  fixed  orthogonal 
aurfaoaa.  Tha  superposition  of  tha  velo¬ 
city  conponents  in  the  above  etreamlike 
function  definitions  gives  tha  three  gas 
valocity  conponents  v^,  Vj  and  V},  which 

identically  satisfy  tha  conservation  of 
aass  equation  (1) 
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In  the  above  aquation*  Qj,  and  Qj  ar* 

three  vorticity  cowpooanto  which  can  be 
expreaaed  in  tarn*  of  the  velocity  deri¬ 
vative*  froe  the  definition  of  the  vorticity 
vector  (equation  •). 

|  the  elliptic  governing  equation*  for 
the  streaalike  function*  xh»  xv  «nd  xc  can 

bn  obtained  free  the  substitution  of 
equations  (2)  through  (5)  into  the  Xj  and  Xj 

ooaponenta  of  equation  ())  and  the  x, 
conponent  of  equation  (•)  respectively. 

!  hq  lifaq  fej  (h3*hn  lKpsj 

-  (13) 

'  hq  IE“KJ  1  *  hq  hqlh2*v] 1 


it  to 


at  gradients 

flow  only  when  it  is 


by  gradients  in  the  velocity 
such  that .gradient  la  stsgna- 

M  aclltaU.  .  In  Mb 


ftion  pressure  exists1*.  Therefore,  in  the 
pcseant  analysis  thn  acxsantuu  equation  is 
sapaeasad  in  taras  of  the  total  pressure, 
I,  Shiah  to  neat  as  ana  of  the  flow  verl- 
•htos  -to  '  the  solution. 


f.n.t,  .  £  jl  ,  E  J» 


■holic  equations  can  be  derived  for 
i  pressure  P  and  the  through  flow 

m  aCuatlon  nr  the  total  pressure 
htalned  through  taking  the  dot 
of  the  velocity  vector  with 
>  CO 


f.7»  ■  0  (10) 

the  following  equation  can  also  be  derived 
for  the  through  flow  vorticity  developaent 
free  the  nensntun  equation* 
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St  is  clear  that  the  above  aquation  also 
represents  the  x,  ooaponent  of  Helnholts 
agnation  uhich  etn  be  written  is  the 
fnllewiag  tern  for  inviseid  coapreesible 
flew 


-  h1h3cT 


fej  IEJKJ  lx^(hlxc,  1  +  hq  lspq  Sxj<hi*c)  1 
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Results  and  Discussion 
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Tbm  analytical  formulation  results  In i 
WMbsmbml  Poisson's  equations  for 
tbs  stress! Its  functions.  Tbsse  equations 
are  Obtained  from  satisfying  the  equations 
of  notion  on  throe  sets  of  orthogonal 
aerfaoee  represented  by  constant  values 
of  the  coordinates  *1#  x2  and  x3  as  shown  ! 

in  Pig.  1.  The  source  tern  in  the  result¬ 
ing  equations  are  dependent  upon  the 
variation  of  the  flow  properties  and  on 
the  flea  in  the  direction  normal  to  these 
surfaces,  bn  iterative  procedure  is  used 
in  the  numerical  solution  because  of  the 
ispsndsncy  of  the  source  ten  in  each  set  1 
of  strosnlihe  function  equations  on  the 
solutions  obtained  for  the  remaining  two 
sets. 


the  iterative  procedure  for  the  numer- 
ioal  computations  consists  of  the  use  of 
a  marching  technique1*  in  the  solution  of 
equations  (7) ,  (10)  and  (11)  along  the 
through  flow  direction  x, ,  and  a  successive 
over  relaxation  method  for  the  solution  of 
the  two  dimensional  elliptic  aquations 
(111  *  (14)  and  (IS)  for  the  streaelike 
functions.  The  flaw  density,  which  is 
allowed  to  lag  one  iteration  in  the  numeri¬ 
cal  computations  is  determined  from  the 
local  total  pressure,  total  stagnation 
enthalpy,  and  from  the  flow  velocity! 
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IT 


(1 


V*  (l/Y-1) 

Iff1 
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The  boundary  conditions  for  the  streaalike 
functions  are  carefully  chosen  to  insure 
the  uniqueness  of  the  solution.16  Dirichlet 
boundary  conditions  are  determined  for  the 
str manlike  functions  from  the  requirement 
of  aero  flux  at  the  duct  boundaries^7  The 
Xj  derivatives  of  the  flow  velocities  are 

set  equal  to  sero  at  the  duct  exit,  while 
the  Oirichlet  boundary  conditions  for  the 
streenlike  functions  at  x,  *  0,  are 
expressed  in  terms  of  integral*  of  the 
inlet  flux.  The  total  pressure,  total 
tanparature  and  through  flow  vorticity 
profiles  are  required  at  the  duct  inlet 
to  start  the  aerching  solution.  Mors 
details  about  ths  numerical  procedure  can 
be  found  in  reference  16. 


The  results  of  the  numerical  compu¬ 
tations  are  presented  in  an  accelerating 
rectangular  elbow  and  compared  with  exist¬ 
ing  experimental  measurements  19 .  The  duct 
was  dssignad  by  Stanits,  0  using  inviscid 
incompressible  two  dimensional  analysis  to 
avoid  boundary  layer  separation.  This 
experimental  data  was  also  used  for  com¬ 
parison  with  results  of  the  numerical 
analysis  by  other  investigators. 15 »21 
The  experimental  maasuraments  wara  obtained 
for  different  inlet  total  prasaura  profiles 
to  investigate  the  effects  of  secondary 
flow.  Tha  inlat  total  prasaura  profilaa 
were  generated  using  perforated  plates  of 
dif farant  haights  as  spoilers.  Figurs  2 
shows  tha  duct  geometry  and  the  computation¬ 
al  grid  in  the  x^  and  x2  directions.  Due 

to  symmetry,  the  flow  computations  were 
performed  in  the  lower  half  of  the  duct 
using  s  (9  x  13  x  S3)  grid  in  the 
x3,  *2  *nd  xi  directions,  respectively. 

Figure  3  shows  the  inlet  velocity  profiles 
which  were  used  in  the  numerical  calculations 
w'ith  no  variation  in  the  xa  direction. 

The  experimental  measurements  for  2.S  inch 
spoillsr  which  were  obtained  half  way  be¬ 
tween  tha  pressure  and  suction  surfaces  ere 
shown  in  tha  same  figure.  The  results  of 
the  computations  ere  presented  in  non- 
dimensional  form  in  Figs.  4  through  8.  The 
flow  voloeitios  arc  normalized  with  respect 
to  the  mexianns  inlet  velocity,  VImh>  while 

the  pressures  ere  normalized  with  respect 
to  the  critical  pressure,  which  corresponds 
to  a  tank  gauga  pressure  of  20  inches 
water  in  the  experimental  measurements. 

The  orthogonal  curvilinear  body 
fitted  coordinates  for  the  flow  computations 
were  generated  numerically  using  the  code 

developed  by  Davie22  which  is  based  on  the  - 

;  Schwarts-Kristoffel  transformation.-  The' 
values  of  the  orthogonal  coordinate  in 
the  transformed  plane  were  between  0  and  1 
in  the  x2  direction  and  0  and  6.791  in  the 
Xi  direction.  The  computational  grid  is 
uniform  in  the  x,  direction  where  half  the 
duct  height  la  equal  to  1.875  times  the 
duct  exit  width.  The  convergence  of  the 
numerical  solution  was  fast  with  CPU  time 
of  3.5  minutes  for  25  outer  iterations 
on  AMDAHL  370.  The  inner  iterations  did 
not  exceed  25  in  all  the  iterative  solu¬ 
tions  for  xh#  Xv  and  xc  solutions  on  all 

77  surfaces. 

The  results  are  presented  first  for 
the  computed  static  pressure  coefficient 
distribution  over  the  duct  curved  walls 
and  are  compared  with  the  experimental 
measurements.  Complete  spanwise  atatic 
prsssura  distributions  over  the  pressure 
and  auction  aurfacas  of  the  elbow  were 
reported  by  Stanits  et  al  19.  for  the  flow 
with  nominal  exit  Mach  number  of  0.26. 

Figures  4  and  5  show  the  computed 
static  pressure  coefficient  distribution 
over  the  duct  curved  boundaries  at  two 
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thedsct  edt.  Similar 
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od  to  via- 
■  is  con- 
ieh  show 

flow  bsglna  to  coll  op  half 
da  duet  length  towards  tba 


oai  ti.  Kao* 
to  m  •,  2.0, 


thus  tbs  oeayutod  total 
are  at  the  cross  sectional 
wondlag  to  «,  -  1,  23,  31 
osasa  sections  eonraspoad 
3.0  and  4.3  in  Stanits 
■  follow  the  development 
f  flow  through  tba  rotation 


while  at  tba  dost 


fftfaao  ?  shews  the  otraaawlsa  vorti- 
sitr  seats ars  at  difforant  cross  sectional 
sartaaas  samspeaitlng  to  a,  •  33,  31 
am  <3  respectively.  the  streaawis* 
vernal  ty  it  sanattisaa  with  respect  to 
m  ratio  of  the  aariaaa  inlet  velocity 
he  tba  dost  exit  width,  the  computed 
rsaalta  pradlat  an  initial  inerease  in  the 
ttaaafh  flow  vortioity  ccaponents  in  the 
mat  half  of  the  dost  length  frost  its 
aaao  veins  at  inlet,  than  a  decrease  in  Ox, 
towards  the  duct  salt.  The  region  of 
■aafchaa  through  flow  vortioity  near  the 
sad  wall  la  soon  to  novo  across  the  pas- 
sege  free  the  pressure  surface  towards  the 
aaaiion  surface,  tt  is  interesting  to 
observe  that  through  flow  vortioity  daval- 
ocMBt  nmAdfl  loroii  tht  pumh  And  it 
aat  Halted  to  the  end  wall?  this  is  not 
aarprlalag  since  one  can  see  froa  rig.  3 
that  tha  shear  flow  occupies  shout  half 
tba  area  at  tha  duet  inlet. 


k  aothod  is  praauatod  for  the  aaalysia 
of  latataal  throe  dimnslonsl  ceaprasaibla 
lasted d  rotational  flow.  The  straaalike 
faaflioa  fOraalatlon  loads  to  a  very' 
aaaaeaioal  alllptie  solution,  that  is 
sal  fade  for  tarboaaehlnary  applications. 


Tha  aaalyais  is  ganeral  and  applicabla  to 
flow  flalda  with  both  total  prsssura  and 
total  taaporaturo  gradianta.  Tha  rasults 
of  the  computations  ars  prasantad  for  tha 
flow  in  an  accelerating  rectangular  sibow 
with  shear  ialot  velocity  profile  and  90* 
turning  angle.  The  analysis  prsdiets  tha 
secondary  flow  developaaat  and  the  computed 
results  ars  in  agrssasnt  with  ths  sxperi- 
aaotal  asssurmants  in  tha  regions  where 
tba  viscous  dissipation  effects  era  not 
aignif leant. 
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PRESSURE  SURFACE 


(1) 


ns 


flow 

to  bo  dependent 

tho  difference  schema  for  tbo  continuity 
s tally  la  tho  prediction  of 
flow  developwant.  Tho 
Lotion  fluttered  on  tho  iaoido 
■,  suffered  o  loss  of  ansa  and 
to  folly  oonwor«o,ll  and  tho  author 
that  ooro  attention  to  tho  non- 
linearities  of  tho  flow  nay  possibly 
alleviate  tho  loot  two  problems.  Other 
nothndo  dove loped  for  tho  olliptie  aolution 
of  the  otoody  ate to  full  Nevisr-Stokes 
or*  baaed  on  tho  ortonaion  of 


tho.1 


tfe*  well-known  3-0  atroaai  fuaotiog  vortieity 

tmmi ia|ion  to  throw  d  logo  along.  *3.14 

mS  and  iolluawlJ  defined  a  vector  poton- 
tiol  to  idontically  oatiafy  the  aquation 
of  oonaezvatidn  of  naaa  in  thro*  dinenaiooal 
flow  field*.  Tho  vector  potential  vorti- 
dty  fomnlation  has  boon  used  to  aolvo 
problem  of  looinor  natural  oonvoe 
t.13,14  |g  this  foraMlatian,  the 

governing  aquations  consist  of 
3-0  Poisson  equations  for  the  voetor 
potential  and  throe  vortieity  transport 
agnations,  no  aain  advantage  of  this 
fomnlation,  in  two  dtawnalona,  nanely  tho 
await or  sis*  of  the  governing  equations 
is  actually  reversed  in  throe  dimensions. 
OtfOe  throo-diaanoional  differential  equa- 
tlons  have  to  bo  solved  for 
toe  thro*  oonpononto  of  tho  vector  poten¬ 
tial. 


On*  oon  conolud*  fro*  the  preceding 
diacussione  that  existing  space-elliptic 
salvers  of  the  3-0  ■avier-dtohes  equations 
am  very  oostly  la  to ms  of  CPU  tine  and 
stoeage  requirenants.  In  addition,  sane 
of  these  netted*  which  were  developed  for 
sinple  convection  problens  have  not  bean 
very  successful  in  through  flow  calcula¬ 
tions.  The  presented  work  represents  a 
new  formulation  for  the  3-0  Havier-Stokes 
equations  that  leads  to  an  alliptic 
solution.  The  fomulation  is  based  on  the 
us*  of  the  2-0  streamlike  functions15  to 
identically  satisfy  the  equation  of  conser¬ 
vation  of  aass  for  3-0  rotational 
flows.1*'17  The  governing  equations  con¬ 
sist  of  the  vortieity  transport  equation 
and  2-0  Poisson  equations  for  the  stream¬ 
ings  functions.  The  present  netted  is 
very  general,  in  that  inviscid . f low  solu¬ 
tions  can  be  obtained  in  the  Unit  when 
hs  -*  «.  Zn  fact,  numerical  solutions  have 
been  obtained  for  inviseid  rotational 
inoonpreoaible1*  and  compressible17  flows 
in  curved  duets  and  it  was  demonstrated 
that  the  method  can  predict  significant 
seoondary  flow  and  straanwise  vortieity 
development  due  to  inlet  vortieity.  The 
following  work  represents  the  generalisa¬ 
tion  of  this  fomulation  to  internal  three- 
dimensional  viscous  flow  problems. 


Analysis 


fora,  for  iaeospressible  viscous  flowt 

<$*y>5  -  -  i-  g2a 


-  0  (2) 

Zn  the  above  equations  a  -  V*D/v  when 
the  velocities  are  normalised  by  V*,  the 
space  dimensions  by  O  and  the  vortieity  by 
V*/D. 


The  solution  to  the  three  dimensional 
viscous  flow  is  obtained  in  tens  of  the 
three  vortieity  oonpononto  and  three  stream- 
like  function*! 5  which  are  defined  to 
identically  satisfy  the  equation  of  con¬ 
servation  of  aass  for  general  three- 
dinansional  rotational  flow  fiolds.  Unlike 
the  traditional  stroaa  function  solutions 
whlah  must  b#  obtsinsd  on  stroaa  surfaces, 
the  following  strsamliks  function  velocity 
relations  permit  the  definition  of  these 
two  dimensional  functions  Xi(*,y),  x2(y**)» 

Xj(x,s)  on  fixed  non-stream  surfaces  in 
the  flow  field. 


ftreamlike  functions  Velocity  halations 

The  streenlika  function  formulation 
was  developed  by  the  authors15  to  model 
Internal  three  dimensional  flow  fields.1* 
More  details  and  general  definitions  in 
curvilinear  coordinates  of  the  strsamliks 
function  can  bo  found  in  reference  17. 

Por  the  sake  of  simplicity,  the  equations 
will  be  presented  here  for  incooprasaible 
flow  using  Cartesian  coordinates. 


Definition  of  x^ 
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The  governing  aquations  oonsist  of  ■  Ju, 

the  vortieity  transport  equation  and  the  .  3  .  w  _  7 

equation  of  conservation  of  aass,  which  *x  3  *  ®x 

are  written  in  the  following  dimensionless  0 


(5b) 
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The  •lliptio  lyittn  of  equations  are 
solved  using  aa  Iterative  procedure.  At 
each  global  1 tar at ion,  tha  linaar  aquations 
were  solved  by  suecassiva  relaxation 
methods.  numerical  computations  in  a 
straight  duet  with  £/DRe  >0.1  vara  par- 
formed  using  a  uniform  grid  with 
Ax/DRa  -  As/DRs  -  0.001  and  Ay/DRs  -  0.0033. 
Oua  to  symmetry ,  tha  computations  wars  only 
carried  out  in  one  quarter  of  tha  duct  for 
X1#  X)»  C  and  {  since  c(x,y,s)  -  -n(s,y,x) 

and  Xj_ <*»y)  *  -x2  (*»y) .  Relaxation  para¬ 
meters  of  1.4  for  x^»  l.f  for  X3  and  0.4 

for  c  and  {  wars  used  in  tha  inner  itera¬ 
tions  with  a  convergence  criteria  of 

«-  ■  1  x  10“*,  c,  ■  1  x  10”5,  cv  -  lxio"5. 


*x  -  3  x  10 
aquation 1 


according  to  tha  following 


C  »  • 


i  1  *  “i#j|>  - 


1  *1** 


9m  Xb«t  of  tbs  boundary  conditions 
das  to  agatHtry  at  the  planes  x  •  0 

At  x  •  0 


I  •  C  •  0 


x2  -  x,  -  • 


At  S  -  0 


a  -  C  -  0 


xa  -  Xj  -  0 

fully  developed  flow  conditions 
applied  at  the  duet  exit. 


C  -  0 


gi.jja., 


The  numerical  solutions  required  50 
global  iterations  and  a  CPU  time  of  2 
mlnmtaa  and  13  seconds  on  AMDAHL  370 
using  an  11  x  11  x  14  uniform  grid. 

The  overall  number  of  iterations  was  357 
far  tbs  vortioity  equations  and  179  for  tha 
stream! lira  function  equations.  Tha  numeri¬ 
cal  solution  dcmisln  extended  1.C7  diameters 
upstream  of  tbs  duct  inlat,  where  tha  flow 
velocity  was  taken  to  ha  uniform  and  equal 
to  one.  The  results  of  the  numerical 
computations  axe  presented  at  y/DAa  -  0.0, 
0.01  and  0.10.  The  through  flow  contours 
at  the  duet  inlet  are  presented  in  Pig.  2. 
The  flow  development  from  a  uniform 
through  velocity  to  the  profile  of  Pig.  2 
at  the  duct  inlet  is  accompanied  by  lateral 
flow  displacements  due  to  the  secondary 
velocities.  The  secondary  velocity  con¬ 
tours  at  the  duct  inlet  are  shown  in 
Pig.  3  for  the  vertical  velocity  component 
w.  The  elllptieity  of  the  numerical  solu¬ 
tion  is  demonstrated  in  the  velocity  con¬ 
tours  at  the  duct  inlet,  and  in  tha  velo¬ 
city  fields  up  to  0.13  diameters  upstream 
of  the  duct  inlet,  figure  4  shows  the 
contours  for  the  secondary  velocity  com¬ 
ponent,  w,  at  y/DRe  “  0.01.  A  comparison 
of  Pigs.  3  and  4  reveals  the  change  in 
both  the  magnitude  and  the  location  of  tha 
■axlmy  secondary  velocities  along  the 
duet.  The  development  of  the  secondary 
velocity  component,  w,  along  the  plane- of 
symmetry,  x  •  0,  is  presented  in  Pig.  5. 

One  can  see  that  the  maximum  secondary 
velocities  are  found  near  the  solid 
boundaries  at  the  duct  inlet.  As  the  flow 
proceeds  towards  fully  developed  conditions, 
the  secondary  velocities  decrease  and  the 
location  of  the  maximum  values  moves 
toward  the  center  of  the  duct.  The  re¬ 
sults  of  ths  numerical  computations  at  the 
duct  exit  are  shown  in  Pig.  6  for  the 
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Fig.  2.  Through  Flow  Velocity  Contours  at 
tha  Duct  Inlat  (y/DRa  •  0.00) . 


Worticity,  c  Contours  at  ttaa  Duet 
Kntranc*  (y/DHs  •  0.0) . 


tig.  9.  Straaalika  Function,  x3  Contours 
at  ths  Duct  Xnlat  (y/DRs  -  0.0). 


tun  using  priatttlfi  rariauiee.  xa  tax* 
UmlA  is  oouputad  fro*  tha  solution  of  a 

_ _ _  with  —min  boundary  conditions.  Both 

Eofr  therefore  nquirt  u  solution  of  throa  para- 
'“‘“'“Hons.  The  first  nsthod  requires  in 
of  ths  3-0  Poisson  aquations  with 


aloof  tha  third  boundary,  whila 
3-0  Poisson  equation  for 
conditions  over  ail  tha 
tha  relative  esrits  of  these 
of  CPO  tine,  oonputar  storage  require- 
tiao  whan  iterative  and  direct 
ths  solution  of  Poisson  aqustion. 

solution  for  viscous  flow  in  a 
~  an*  curved  duct  (101  frost  the  governing 

variables  using  a  new  finite  diffaraaea 
Sira  encountered  in  tha  application  of 
f| thveogh  How  calculations;  tha  tandaney  of  tha 
00  flutter  and '  a  loss  of  ansa  bateaaa  tha  duct  inlet 
issulfes.  .Nora  recently.  Dodge  (11] 

^  Sitoiiiy^Oplitt-iKOoadura  in  which  tha  velocity 
vleocue  and  potential  conponants,  and  tha 
field  is  dotsmlnad  front  tha  potential  velocity 

aquation  for  the  viscous  velocity 
ia  obtained  Croat  tha  ssnantun  aquation  with 
in  tarns  of  tha  derivatives 
.ty  potential  vector  oonponeat,  while  the  govern* 
for  the  velocity  potential  vector  component 
continuity.  Beyond  thin  formulation, 
solution  procedure  is  partially  parabolic 
h*  neglected  tha  streaawise  diffusion  of  momantw  In 
a  narabing  solution  for  the  viscous  velocity  coat- 
analysis  Itself,  in  tares  of  the  type  of  the 
felons  sad  their  boundary  conditions,  is  com¬ 
parable  to  thn  velocity  pressure  formulation  since  tha  govern¬ 
ing  equation  for  the  velocity  potential  is  a  3-D  Poisson 
aquation.  Dodge  did  not  discuss  the  boundary  conditions  for 
tba  velocity  potential.  Be  only  aantioned  that  it  can  be 
coqplem  end  that  he  used  aero  potential  gradient  normal  to 
the  wall  ia  his  inam’ioal  solution,  other  formulations  that 
osa  land  to  elliptic  solution,  wars  described  in  references 
(SI  sad  (•] ,  however  they  will  not  be  discussed  here  since 
they  have  not  yet  been  applied  to  3-D  flow  computations. 


Xa  sunnary,  existing  elliptic  solvers  for  the  Bavier- 
•tokes  equations  require  the  solution  of  one  or  thraa  3-D 
Poisson  equations  in  addition  to  tha  thraa  momentum  or  tha 
three  vortieity  transport  equations.  In  the  case  of  the 
velocity  pressure  formulation,  the  velocity  vector  ie  evalu¬ 
ated  from  the  momentum  equation  and  tha  continuity  aquation 
is  satisfied  indirectly  in  the  pressure  equation.  The  con¬ 
vergence  Of  the  iterative  numerical  procedure-  would  be  very 
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tea  IHywlil  atetoar  (M  •  ,  appaara  In  tha  •qottlooa  as 

I  wwu  of  iwmaHttng  tte  valoeitias  by  v*,  tte  ipto« 
iteatelans  by  P  and  th*  oortieity  by  v*/D,  wtera  v  is  tte 


act  written  in  Cartaaian  coordlnatas  aa 


•ilpflipi* 
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afeava  n,t,  C  an  tea  ooaponanta  of  tha  vortieity  vaetor,  S,  in 
R|*  S*  y  and  a  directions,  respectively;  u,v,w  ara  tte 
awapMMMHi  Of  tte  velocity  factor  V. 


fertieity  Velocity  Kalations 
11  *  fy  "  la 

0 

c  .3 s.ia 

(  is  >r 


three  atgaaalite  functions  [12]  xh(x,y) ,  ^(y**)  #  Xc<x#*) 

an  dafiaad  to  idantieally  satisfy  tha  aquation  of  conserva¬ 
tion  of  aaao  in  tea  eaaa  of  intarnal  rotational  flows  [13, 

14].  the  velocity  fiald  is  datarainad  from  tha  streaolike 
functions  according  to  tte  following  rslations i 


riym  3  Hww  tbs  boundary  conditions  for  the  straw 
lika  fnmtlnns  at  the  extension  of  ths  duet  boundaries, 
tbs  conditions  st  tbs  planes  of  syssotry  arc  unchanged. 


Figure  3 


relations  satisfy  tha  ialst  condition 
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flow  conditions  arc 
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Figure  5 


3.  mmuct  MB  DXSC08SX0H 

The  results  of  the  m— ricel  computations  ere  presented 
in  one  quarter  of  a  square  duct.  Xn  this  case,  only  two  of 
the  vorticity  equations  teqs.  (4)  and  (5) ]  and  two  of  the 
streanlike  function  equations  teqs.  (19)  and  (20)]  are  solved, 
sines  Xjj  *  *ad  C(x#y,a)  ■“  n(s,y,x) .  Referring  to  Fig.  1,  the 

solution  was  obtained  for  Re  ■  50  in  a  duet  with  VDRe  »  0.1 
and  La/DRe  *  0.01  using  30R  and  a  (11  x  11  x  34)  grid. 

Figures  6a  and  6b  show  the  through  *flow  velocity  contours 
at  the  duet  entrance  and  exit.  The  influence  of  the  cascade 
entering  on  the  elliptic  solution  is  demonstrated  in  contours 
of  Fig.  6a. 

lbs  contours  fbr  the  secondary  velocity  couponent,  w, 
are  shown  at  y/DRa  ■  0.0  and  0.0075'  in  Figs.  7a  and  7b. 

From  these  figures,  one  can  see  a  large  change  in  both  the 
sagnitede  and  the  location  of  the  aarletmi  secondary  velocities 
along  the  dust.  The  development  of  the  through  flow  velocity 
profiles  along  the  plane,  of.  symmetry ,_.x  «_-0,_is_shavn.in. 


c-fft. 


rif.  1  for  the  counted  nnltv  and  th«  experimental  amasure- 
manta  ©£  rtftnoct  (16] .  Om  can  that  tha  computed 
Units  act  in  good  agreement  with  tha  experimental  measure- 
agnte,  at  y/Wle  •  0.0075  and  0.02,  hot  that  tha  computed 
"nifiatfh  'flow  velocities  at  y/I»a  -  0.1  did  not  raaeh  tha 
measured  tally  developed  profile. 


DATA  (Hat.  (16]). 


fly.  I.  DmlofMot 


at  tha 


ant  of  tha  through  How  Velocity  Profile 
Central  Plana  x  -  0. 


figure  9  shows  tha  developnent  of  tha  secondary  velocity 
aaaponaot,  w,  along  tha  duet  plana  of  symmetry,  x  -  0. 

Me  experimental  aeasureneats  are  available  for  conparison 
with  tha  ccnputad  secondary  velocities.  Figure  9  shows  that 
the  maximnm  secondary  flow  is  initially  located  near  the 
solid  boundaries,  than  novas  towards  tha  canter  of  tha  duct 
and  decreases  as  tha  flow  proceeds  towards  fully  davalopad 
conditions,  tha  ccnputad  through  flow  velocity  development 
along  the  duet  centerline  are  compared  with  tha  experimental 
measurements  of  reference  [16]  in  Fig.  10.  One  can  sea  in 
this  figure  that  tha  elliptic  solution  predicts  an  increase 
in  the  centerline  through  flow  velocity  in  the  cascade  entry 
region  preceding  the  actual  duct  entry.  Figure  10  shows  that 
tile  computations  slightly  underestimate  the  centerline 
velocity.  Considering  the  coarse  grid  used  in  the  numerical 
computations  ( (11  x  11  x  34)  grid  points  in  the  duct] ,  the 
agreement  of  the  coaputad  results  with  the  experimental  data 
as  shown  in  Figures  8  and  10  is  very  satisfactory. 
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rig.  9.  The  Developnant  of  Secondary  Plow  Valoclty  Profile 
at  the  Central  Plana  x  -  0.  , 


0  EXPERIMENTAL  DATA  (REF.  [16]) 
2.2r  — -  PRESENT  RESULTS 


presents  a  fast  efficient  method  for  tho  3-0 

Barier-Stoke*  aquation*,  it  is  baaad. 
ha  funotlon  rortlcity  formulation  which  toads 
equations  with  Oirichlot  boundary  conditions 
typist  Ur, s  functions,  in  addition  to  the 
past  aquations.  Tha  aathod  is  nora  aconoadcal 
PJtoWr.ftohas  elliptic  solvers,  yat  does  not 
IKiM  of  the  parabolisad  and  paratiatly  para- 
Ms*  It  offers  a  useful  tool  for  tha  nnaarical 
iatarnal  viscous  flow  fields  where  surface 
awastraaar  affects  are  significant,  as  in  turbine 
>t  :Sha  results  of  the  computation*  are  presented 
-&m  in  a  constant  area  duet,  as  a  corner 
panersl  future  applications. 
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Cincinnati,  Ohio  45221 


»  waits  of  as  super  1—ntal  investi- 
«C  tfca  three  <j— m atonal  flew  devel- 
.  fat  •  highly  carved  daet  with  inlet 
profile  ere  presented.  The  three 
.. lasts  of  the  sir  velocity  in  a  curved 
with  e  root angular  cross  section  are 

- jsed  using  Laser  Doppler  Anaaosatry. 

•ignifieaat  through  velocity  oontour  rota* 

— - 1  reported  with  secondary  velocity 

sgnltudas  up  to  0.25  of 
ity  in  the  90*  turning 
dnot. 


The  eacaSdary  flow  affects  the  overall 
tailhaaeahlaary  perforsance  through  its 
influance  on  the  angle  and  the  energy  die* 
trihetien  of  the  flow  leaving  the  blade 
case*  A  deviation  in  thn  exiting  flow 
angina  froa  these  predicted  by  the  blade 
alasaet  aaalysia  results  froa  the  secondary 
uadoeitlae.  The  secondary  losses  are  also 
CSeasd  by  the  redistribution  of  tbe  low 
saergy  flow  by  ths  aea  secondary  veloci¬ 
ties.  These  arose  velocities  ere  associ¬ 
ated  with  tbe  secondary  velocity  develop- 
sent  in  the  a  tree— i—  direetion  through 
the  turning  of  the  flow  with  nonuniform 
inlet  conditions  in  the  blade  paesagea. 

Secondary  flow  in  eoapraaaor  and 
turbine  oasoades  baa  been  the  subject  of 
several  theoretical  and  experimental 
investigations.  In  most  of  the  experi¬ 
mental  secondary  flow  investigations,  the 
flew  aaeoaraaoete  hove  been  limited  to 
cascade  inlet  and  exit  conditions  to  pro¬ 
vide  anpirieal  correlations  for  secondary 
flew  loenee  and  exit  flow  angles. 

Laagaatan*  obtained  detailed  neesure- 
aaata  showing  the  general  characteristics 
ef  the  end  wall  flew  downstream  of  the 
saddle  point  and  upstream  of  the  turbine 
eeacade  passages.  Moore2  measured  the 
•sent  pressure  and  flow  directions  down¬ 
stream  of  the  turbine  cascade  trailing 
edge.  Oeveral  investigators  measured  the 
development. of  secondary  flow  la  curved 

duets4*7  from  e  fully  developed  inlet 


•  Professor;  Associate  yellow  axaa. ' 
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velocity  profile.  These  studies  demon¬ 
strate  the  seeondary  flow  developawnt  in 
the  absence  of  added  complexities  of  the 
cascade  blade  leading  edges. 

In  the  present  work  detailed  — asure- 
aant  of  tba  three  flow  velocity  components 
are  obtained  in  e  curved  duct  with  a  nearly 
linear  shear  flow  inlat  profile  produced 
using  a  grid  of  parallel  rods  with  varying 
■pacing.  Under  these  conditions,  the 
development  of  the  secondary  velocities 
associated  with  the  passage  vortex  is  not 
lied  ted  to  the  boundary  layer  region  near 
the  wall,  but  extends  instead  through  the 
whole  pea saga  sections.  The  expert— ntal 
measurements  of  this  complex  flow  field 
are  baaed  on  the  use  of  a  two-color  back 
scatter  Laser  Doppler  Veloeimater. 


EXPERIMENTAL  SET-UP 

Tbe  a xperi— ntal  set-up  la  shown 
■ehaaatlcally  in  Pig.  1.  It  consists  of 
the  tunnel,  the  seeding  particle  atomiser, 
the  LDP,  optical  and  data  acquisition 

systems. 

Tunnel 

The  high  pressure  air  supply  from  the 
storage  tanks  ia  regulated  to  a  lower 
pressure  before  entering  the  12"  diameter 
settling  chamber.  A  1.5"  thick  honeycomb  of 
0;187"  cell  and  0.003"  wall  thickness  is 

placed  in  a  4"  die— ter  p VC  tube  to  , - 

condition  the  flow.  The  latter  extends ' 

18”  Inside  the  chamber  and  blsnds  smoothly 
into  a  22.75"  long  rectangular  channel 
preceding  the  curved  duct.  The  duct  is 
shown  schematically  in  Fig.  2  and  consists 
of  a  90*  band  of  6"  —an  radius  and  a 
2"x4"  rectangular  cross-section.  The 
duct  walla  era  made  of  plexiglass.  The 
thickneaa  of  the  curved  wall  ia  equal  to 
1/8"  while  tba  plane  walla  are  3/4"  thick. 

The  curved  duct  ia  connected  to  straight 
ducts  downstream  and  upstream  where  the 
sheer  flow  ia  produced  using  a  grid  of 
parallel  rods  with  varying  spacing.  The 
grid  imposes  a  resistance  to  the  flow  that 
varies  across  the  section  ao  as  to  produce 
variation  in  the  flow  velocity,  without 
introducing  an  appreciable  gradient  in 
atatlc  pressure,  me  basic  relations 
between  rod  also,  spacing  and  the  produced 
velocity  gradient  were  derived  first  by 
Owen  and  Zienkiewicz* »and  modified  later 
by  Llvmaay  and  Turner9  and  by  Elder10  for 
more  generalised  profiles.  The  basic 
relations  for  a  grid  to  generate  a  uniform 
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evrvad  wall.  Finally  tha  mountings  of  the 
sanding  and  raoaiving  optios  axa  dasignad 
to  allow  for  rotation  around  tha  optical 
aria  op  to  t0*»  in  ordar  to  obtain 
tha  aaaauraaianta  of  tha  two  valooity 
components  in  any  apaeifiad  dir action. 
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An  signal  processors  and  a  digital  Mine 
11/83  eonpatere  wars  used  on  line  to  acqu¬ 
ire  synchronised  data  for  tha  siaultaneous 
notour awanta  of  the  taro  velocity  components 
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The  high  prassurs  air  supply  fron  air 
storage  tanks  was  regulated  to  a  pressure 
Of  2.0$ 84  x  10*  H/a2  gauge  (30  psig)  at  the 
orifice  nater,  using  a  pressure  regulating 
valve,  to  give  an  air  nass  flow  of  0.1143- 
kg/sae  (0 . 252  lb/sac) .  This  corresponds 
to  Xeynolda  nunbar  of  1.3  x  105,  based  on 
the  height  of  the  duet  and  the  nean  inlet 
velocity  of  20  a/ sec  ($6  ft/sec) ,  upstream 
of  tha  shear  velocity  generator. 

Tha  axperinental  velocity  measurements 
were  obtained  at  sections  B,  C,  D  and  s  as 
shown  in  Fig.  2.  Sections  C  and  D  are 
located  in  tha  curved  duct  at  the  45°  and 
75°  turning  angles,  while  the  first  and 
last  naaauring  stations  B  and  8  are  located 
in  the  straight  portions  of  the  tunnel.  A 
quarter  inch  spacing  between  the  measuring 
points  in  tha  radial  direction  and  also  in 
tha  direction  nornal  to  tha  duct  plane 
walls  was  kept  in  tha  neasuranants  at  all 
four  cross  sections.  The  velocity  measure- 
neats  were,  therefore,  obtained  at  105 
points  of  a  7  x  15  grid  in  every  section. 
In  ordar  to  detaraina  the  three  velocity 
oonpooents  at  each  measuring  point  the 
neasuranants  ware  obtained  once  with  tha 
Laser-optics  axis  perpendicular  to  the  duct 
plane  wall,  than  repeated  with  tha  Laser- 
optics  axis  perpendicular  to  the  curved 
wall.  Tha  first  sat  of  measurements 
provided  tha  through  flow  velocity  O,  in 
tha  direction  normal  to  the  tunnel  cross- 
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Motions,  and  ths  radial  velocity  ooapo* 
nants,  v_,  while  the  saooad  sat  of  measure¬ 
ments  fin  tba  third  valoolty  oompoaent  Vy 
and .also  tba  through  flow  valoolty  0.  Tba 
through  flow  velocities  frca  tba  two  sats 
'  of  measurements  wars  oomparsd  to  datarmioa 
tba  rapaatablllty  of  tba  data  aftar  tba 
doct  la  rotatad  90°  relative  to  tba  settl- 
lag  obawbar.-  Tba  dlffaranoa  batwaaa  tbasa 
two  vs  1ms  was  not  found  to  aaoaad  3. St  in 
tba  reported  aaasnrananta.  Tba  raaulta  of 
tba  experimental  measurements  of  tba  tbrM 
valoolty  components  aoraallsad  with  re¬ 
spect  to  tba  naan  flow  valoolty  Ok,  up* 
atraan  of  tba  shear  generating  grid,  ara 
praaaatad  la  Figaros  4  through  11. 


tba  profllas  of  tba  nornallsad  through* 
flow  valoolty  (0/0^)  ara  praaaatad  at  savan  . 
conoantrio  cyolindrical  aurfaoaa  batwaan 
tba  duct  lunar  and  outer  curved  walla  at 
aectiona  B,  c,  0  aad  8  la  flguraa  4,  5,  6 
aad  7  respectively.  Oaa  can  saa  that  the 
abaar  generating  grid  produoad  tba  daairod  • 
lidear  velocity  variation  la  tba  a  direct¬ 
ion  except  la  tba  region  near  tba  upper 
wall  vbera  tba  valoolty  gradient  Is  higher. 
This  deviation  was  due  to  the  influence  of 
tba  diataaca  of  tba  last  grid  wire  frca  tba 
upper  wall.  This  factor  was  not  found  to 
have  significant  affect  near  tba  lower  wall 
where  the  wire  grid  spacings  ara  saaller. 
Other  abaar  valoolty  generating  grids  with 
different  wire  disasters  (0.039*  and  0.125*) 
have  been  investigated.  The  foroaentioned 
effect  wm  even  snra  pronounced  la  the 
case  of  larger  wire  disaster  grid.  On  tba 
other  hand,  tba  velocity  profiels  produced 
by  tba  grids  of  tba  saallar  wire  disaster 
vara  found  to  produce  velocity  variations 
along  the  wire  length.  Careful  examina¬ 
tion  of  tba  grid  revealed  non-uniformities 
in  tba  wire  spacings  la  this  direction, 
which  was  found  difficult  to  control. 
Figures  5  through  7  daaonatrata  the  change  >• 
it  tba  through  flow  velocity  profile  with 
the  duet  turning  angle,  initially,  the  flow 
accelerates  along  the  inner  wall  and  decel¬ 
erates  along  tba  outer  wall  to  approach  po¬ 
tential  free  vortex  velocity  distribution. 
This  can  be  aeon  by  comparing  the  velocity 
profiles  in  Figures  4  and  5,  at  x  ■  1.75” 
and  x  ■  0.25”  raapactivaly.  Later  on  the 
flow  decelerates  along  the  inner  wall  and 
accelerates  along  the  outer  wall.  This 
pattern  is  reinforced  by  the  secondary  flow 
velocity  davelopawnt,  which  tend  to  trans¬ 
fer  the  slower  moving  flow  towards  the  in¬ 
ner  wall  aa  can  be  seen  by  comparing  the 
profiles  at  x  •  1.75”  in  Figures  6  and  7. 


The  pt 

ary  velocity  in  the  radial  direction  (Vx/0h) 
are  also  presented  at  seven  concentric 
cylindrical  surfaces  which  are  equally 
spaced  at  0.25”  between  the  inner  and  outer 
curved  wells  of  the  duct.  Figure  I  combines 
all  the  profiles  at  sections  C,D  and  S,  as 
the  LOV  measurements  did  not  indicate  any 
significant  secondary  velocities  at  section 
8.  On  the  other  hand,  the  profiles  for  the 
normalised  secondary  velocity  for  the 
normalised  secondary  velocity  in  ths  verti¬ 
cal  direction  (V  /t^)  are  presented  at  15 


parallel  planes  which  are  equally  spaced  at 
0.25"  between  the  duct  plane  walls  as  shown 
in  Fig.  9.  The  same  symbols  which  were 
used  for  the  through  flow  velocity  profiles 
at  sections  C,  D  and  8  are  also  maintained 
in  presenting  the  measured  secondary  veloc¬ 
ities  in  Figures  0  and  9.  The  vertical  com¬ 
ponent  of  the  secondary  velocity,  V  ,  was 
not  Masured  at  section  0-0,  due  to  the 
deterioration  of  the  quality  of  the  duet 
curved  outer  wall  after  obtaining  the 
measurements  at  the  other  sections.  The 
maximum  secondary  velocity  components  were 
measured  at  section  B-B  and  were  found  to 
be  0.264  in  the  radially  inward^ direction 
near  the  duct  upper  wall  and  0.138  in  the 
vertically  upward  direction  near  the  duct 
inner  wall.  These  values  are  also  nor¬ 
malised  by  the  inlet  flow  velocity,  ob, 

'before  the  shear  generating  grid. 

The  two  measured  secondary  velocity 
components  were  combined  to  produce  the 
secondary  flow  patterns  at  sections  C  and 
8  which  are  shown  in  Figures  10  and  11, 
respectively.  From  these  figures  the 
developswnt  of  the  passage  vortex  due  to 
the  nearly  linear  shear  inlet  flow  profile 
can  be  observed  throughout  the  duct  cross- 
sections.  The  secondary  velocities  assoc¬ 
iated  with  this  vortex  tend  to  move  the 
slower  flow  in  the  lower  duct  sections 
towards  the  inner  curved  wall  and  the 
faster  flow  in  the  upper  duct  sections,  to¬ 
wards  the  outer  curved  wall.  One  can  see 
from  Oig.  10  that,  the  center  of  the  pass¬ 
age  vortex  at  the  45°  turning  angle  is 
nearly  in  the  middle  between  the  inner  and 
outer  curved  walls  but  closer  to  the  duct 
upper  wall  at  »  -  0.62h.  A  close  examin¬ 
ation  of  Figures  10  and  11  reveals  that  the 
center  of  the  passage  vortex  moves  towards 
the  duct  lower  wall  and  also  in  the  outward 
radial  direction  as  the  flow  turning  angle  ~ 
increases  from  45°  to  90°  between  sections 
C  and  8.  < 
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LDV  measurements  were  presented  for 
the  three  velocity  components  of  the  flow 
in  a  rectangular  curved  duct  with  shear 
inlet  velocity  profile.  Secondary  veloci¬ 
ties  of  magnitudes  greater  than  25%  of  the 
main  velocity,  wars  measured  after  the  90° 
flow  turning  angle.  The  results  demonstr¬ 
ate  tbe  passage  vortex  development  through¬ 
out  the  duet  cross-sections  with  the  flow 
turning  angles.  These  experimental  results, 
oaa  therefore  be  used  to  validate  both 
viscous  and  inviaeid  codas  for  internal 
i  tbrM  dimensional  rotational  flow  fields. 
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